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Two-loop divergences of scattering amplitudes with massive partons 
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We complete the study of two-loop infrared singularities of scattering amplitudes with an arbitrary 
number of massive and massless partons in non-abelian gauge theories. To this end, we calculate the 
universal functions F\ and /2, which completely specify the structure of three-parton correlations 
in the soft anomalous-dimension matrix, at two-loop order in closed analytic form. Both functions 
are found to be suppressed like 0(m 4 /s 2 ) in the limit of small parton masses, in accordance with 
mass factorization theorems proposed in the literature. On the other hand, they are unsuppressed 
and diverge logarithmically near the threshold for pair production of two heavy particles. As an 
application, we calculate the two-loop anomalous-dimension matrix for qq ti near threshold and 
show that it is not diagonal in the s-channel singlet-octet basis. 



I. INTRODUCTION 

Recently, much progress has been achieved in the un- 
derstanding of the infrared (IR) singularities of mass- 
less scattering amplitudes in non-abelian gauge theories. 
While factorization proofs guarantee the absence of IR 
divergences in inclusive observables [1], in many cases 
large Sudakov logarithms remain after this cancellation. 
A detailed control over the structure of IR poles in the 
virtual corrections to scattering amplitudes is a prerequi- 
site for the resummation of these logarithms beyond the 
leading order [2-5]. Catani was the first to predict the 
singularities of two-loop scattering amplitudes apart from 
the 1/e pole term [6], whose general form was only un- 
derstood much later in [7-9]. In recent work [10], it was 
shown that the IR singularities of on-shell amplitudes in 
massless QCD can be derived from the ultraviolet (UV) 
poles of operator matrix elements in soft-collinear effec- 
tive theory (SCET). They can be subtracted by means of 
a multiplicative renormalization factor, whose structure 
is constrained by the renormalization group. It was pro- 
posed in this paper that the simplicity of the correspond- 
ing anomalous-dimension matrix holds not only at one- 
and two-loop order, but may in fact be an exact result of 
perturbation theory. This possibility was raised indepen- 
dently in [11]. Detailed theoretical arguments supporting 
this conjecture were presented in [12], where constraints 
derived from soft-collinear factorization, the non-abelian 
exponentiation theorem, and the behavior of scattering 
amplitudes in two-parton collinear limits were studied. 

It is relevant for many physical applications to gener- 
alize these results to the case of massive partons. The 
IR singularities of one-loop amplitudes containing mas- 
sive partons were obtained some time ago in [13], but 
until very recently little was known about higher-loop 
results. In the limit where the parton masses are small 
compared with the typical momentum transfer among 
the partons, mass logarithms can be predicted based on 
collinear factorization theorems [14, 15]. This allows one 
to obtain massive amplitudes from massless ones with a 
minimal amount of calculational effort. A major step to- 
ward solving the problem of finding the IR divergences 
of generic two-loop scattering processes with both mas- 



sive and massless partons has been taken in [16, 17]. One 
finds that the simplicity of the anomalous-dimension ma- 
trix observed in the massless case no longer persists in 
the presence of massive partons. Important constraints 
from soft-collinear factorization and two-parton collinear 
limits are lost, and only the non-abelian exponentia- 
tion theorem restricts the allowed color structures in the 
anomalous-dimension matrix. At two-loop order, two dif- 
ferent types of three-parton color and momentum cor- 
relations appear, whose effects can be parameterized in 
terms of two universal, process-independent functions F\ 
and /2 [17]. Apart from some symmetry properties, the 
precise form of these functions was left unspecified. In 
this Letter we calculate these functions at two-loop order 
and study their properties in some detail. 



II. ANOMALOUS-DIMENSION MATRIX 

We denote by |.A4 n (e, {p}, {m})) a UV-renormalized, 
on-shell n-parton scattering amplitude with IR singu- 
larities regularized in d = 4 — 2e dimensions. Here 
{p} = {pi,... ,p n } and {m} = {mi, . . . , m„} denote 
the momenta and masses of the external partons. The 
amplitude is a function of the Lorentz invariants Sij = 
Pi • Pj + iO and p\ = mf , where the sign factor 
cry = +1 if the momenta pi and pj are both incom- 
ing or outgoing, and <7jj = — 1 otherwise. For massive 
partons we define 4- velocities Vi = Pi/rrii with vf = 1 
and v® > 1. We further define the recoil variables 
wij = —<Jij Vi ■ Vj — iO. We use the color-space formal- 
ism [18], in which n-particle amplitudes are treated as 
n-dimensional vectors in color space. is the color gen- 
erator associated with the z-th parton and acts as a ma- 
trix on its color index. The product Tj • Tj = T? Tf 
is summed over a. Generators associated with different 
particles commute, and Tf = Cj is given in terms of 
the eigenvalue of the quadratic Casimir operator of the 



corresponding color representation, i.e., C„ = C„ 



for quarks and C g = Ca for gluons. Below, we will la- 
bel massive partons with capital indices (I, J, ... ) and 
massless ones with lower-case indices . . . ). 

It was shown in [10, 12, 17] that the IR poles of such 
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amplitudes can be removed by a multiplicative renormal- 
ization factor Z~ 1 (e, {p}, {m}, A 4 ); which acts as a matrix 
on the color indices of the partons. More precisely, the 
product Z~ X \M.„) is finite for e — > after the coupling 
constant cv^ CD used in the calculation of the scattering 
amplitude is properly matched onto the coupling a s in 
the effective theory, in which the heavy partons are inte- 
grated out [17]. The relation 

Z(e,{p},{m} >t i) = -T({p},{m}, f i) (1) 



dln/x 

links the renormalization factor to a universal anomalous- 
dimension matrix T, which governs the scale dependence 
of effective-theory operators built out of collinear SCET 
fields for the massless partons and soft heavy-quark ef- 
fective theory fields for the massive ones. For the case 
of massless partons, the anomalous dimension has been 
calculated at two-loop order in [8, 9] and was found to 
contain only two-parton color-dipole correlations. It has 
recently been conjectured that this result may hold to 
all orders of perturbation theory [10-12]. On the other 
hand, when massive partons are involved in the scattering 
process, then starting at two-loop order correlations in- 
volving more than two partons appear [16]. At two- loop 
order, the general structure of the anomalous-dimension 
matrix is [17] 
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The one- and two-parton terms depicted in the first three 
lines start at one-loop order, while the three-parton terms 
in the last two lines appear at (D(a 2 s ). The notation 
(j, j, . . . ) etc. means that the corresponding sum extends 
over tuples of distinct parton indices. The cusp angles 
Pi j are defined via 



cosh Pi j 



-sij 



2mimj 



wij . 



(3) 



They are associated with the hyperbolic angles formed 
by the time- like Wilson lines of two heavy partons. The 
physically allowed values are wu > 1 (one parton in- 
coming and one outgoing), corresponding to Pu > 0, or 
wu < — 1 (both partons incoming or outgoing), corre- 
sponding to Pi j — iir — b with real b > 0. These possi- 
bilities correspond to space-like and time-like kinematics, 
respectively. Since in a three-parton configuration there 




FIG. 1: Two-loop Feynman graphs (top row) and one-loop 
counterterm diagrams (bottom row) contributing to the two- 
loop coefficient of the renormalization factor Z s ■ 



is always at least one pair of partons either incoming or 
outgoing, at least one of the wu or vi ■ pk variables must 
be time-like, and hence the functions F\ and ]i have 
non-zero imaginary parts. 

The anomalous-dimension coefficients 7 CU sp (evs) and 
7*(a s ) (for i = q,g) in (2) have been determined to three- 
loop order in [12] by considering the case of the massless 
quark and gluon form factors. Similarly, the coefficients 
7 / (a s ) for massive quarks and color-octet partons such as 
gluinos have been extracted at two-loop order in [17] by 
analyzing the anomalous dimension of heavy-light cur- 
rents in SCET. In addition, the velocity-dependent func- 
tion 7cusp(/3, Oi s ) has been derived from the known two- 
loop anomalous dimension of a current composed of two 
heavy quarks moving at different velocity [19, 20]. 

Here we complete the calculation of the two-loop 
anomalous-dimension matrix by deriving closed analytic 
expressions for the universal functions F\ and /a, which 
parameterize the three-parton correlations in (2). 



III. CALCULATION OF Fi AND f 2 

To calculate the function Fi we compute the two- 
loop vacuum matrix element of the operator O s = 
S Vl S V2 S V3 , which consists of three soft Wilson lines 
along the directions of the velocities of three massive par- 
tons, without imposing color conservation. The anoma- 
lous dimension of this operator contains a three-parton 
term given by 6if abc Tf T 2 b T 3 C Fx(/3 12 , p 23 , p 31 ). The 
function F\ follows from the coefficient of the 1/e pole 
in the bare matrix element of O s . We will then obtain 
ji from a limiting procedure. 

The operator O s is renormalized multiplicatively, so 
that O s Z s is UV finite, where Z s is linked to the anoma- 
lous dimension in the same way as shown in (1). In order 
to calculate the two-loop Z s factor, we have evaluated 
the two-loop non-planar and planar graphs shown in the 
first row of Figure 1, as well as the one- loop countert- 
erm diagrams depicted in the second row. Contrary to 
a statement made in [16], we find that F\ receives con- 
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tributions from all five diagrams, not just from the non- 
planar graph. The most challenging technical aspect of 
the analysis is the calculation of the diagram involving 
the triple-gluon vertex. We have computed this diagram 
using a Merlin-Barnes representation and checked the an- 
swer numerically using sector decomposition [21]. We 
have also checked that for Euclidean velocities our re- 
sult for the triple-gluon diagram agrees numerically with 
a position-space based integral representation derived in 
[16]. Combining all contributions, we find 



Wl2,0!3,/%l) 



12tt 2 ^ 
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where we have introduced the functions 



r{(3) = (3 coth/3, 

/3 2 + 2/3 ln(l-r 
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The function f 2 can be obtained from the above result 
by writing w 23 = -^23 v 2 -ps/ms, W31 = -031 v\ -P3/ m 3 
and taking the limit m 3 — > at fixed Vi-ps. In that way, 
we obtain 
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(6) 

Whether a factorization of the three-parton terms into 
two functions depending on only a single cusp angle per- 
sists at higher orders in a s is an open question. 

It is interesting to expand the two functions r(/3) and 
g{(3) about the threshold point (3 = in — b with b — > + . 
We find 



ITT 



r(/J) = -- + 1 + 0(6), 
9(P) = 



b 
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Based on the symmetry properties of Fi and /2, it was 
concluded in [16, 17] that these functions vanish when- 
ever two of the velocities of the massive partons coin- 
cide. Indeed, this seems to be an obvious consequence 
of the fact that F\ is totally anti-symmetric in its ar- 
guments, while / 2 is odd in its second argument. This 
reasoning implicitly assumes that the limit of equal ve- 
locities is non-singular, but is invalidated by the presence 
of Coulomb singularities in r(/3) and g{(3) near threshold. 
Consider the limit where two massive partons 1 and 2 are 
produced near threshold, with a small relative 3-velocity 
V12 = v\ — v 2 defined in their rest frame. It is then 
straightforward to derive that 

2 

lim h = t~% W 2 + 2*7rln(2|t?i2|)] cos0, (8) 



where 6 is the scattering angle formed by the 3-momenta 
of particles 1 and 3. A similar formula holds for F\. This 
result is anti-symmetric in the parton indices 1 and 2 as 
required, but it does not vanish at threshold. 

Another interesting limit is that of large recoil, where 
all the scalar products wjj become large in magnitude. 
In that case, both F\ and / 2 are suppressed like 0(l/w 2 ), 
because for large (3 
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Note that the non-planar contribution from the first 
graph in Figure 1, which was studied in the Euclidean 
region in [16], contains the leading-power term 



non— planar 
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and is unsuppressed in this limit. However, this contri- 
bution cancels against a leading-power term in the planar 
and counterterm contributions. 

Using that Wu — —su/(2mimj), we see that the 
large-recoil limit corresponds to mimj — > at fixed 
sjj. It follows that the three-parton correlation terms 
described by F\ and f 2 vanish like (vajmj /sjj) 2 in the 
small-mass limit. This observation is in accordance with 
a factorization theorem proposed in [14, 15], which states 
that massive amplitudes in the small-mass limit can be 
obtained from massless ones by a simple rescaling pre- 
scription for the massive external legs. 



IV. ANOMALOUS DIMENSION FOR qq 
NEAR THRESHOLD 



tt 



As a sample application, we apply our formalism to the 
calculation of the two-loop anomalous-dimension matri- 
ces for top-quark pair production near threshold in the 
qq — > tt channel. This matrix (along with the correspond- 
ing one in the gg — ► tt channel) forms the basis for soft- 
gluon resummation at the next-to-next-to-leading loga- 
rithmic (NNLL) order. We adopt the s-channcl singlet- 
octet basis, in which the pair is either in a color-singlct 
or color-octet state. For the quark-antiquark annihilation 
process qi(p\) + q~k{Pi) — ► kipz) + tj{Pi), wc thus choose 
the independent color structures as c\ — Sij Ski and c 2 = 
(t a )ij (t a )ki- In the threshold limit s = 2pi ■ p 2 — ► 4m 2 
it is convenient to define the quantity (3 t — y/l — 4m 2 / s, 
which is related to the relative 3-velocity v t i between the 
top-quark pair in the center-of-mass frame by \v t i\ = 2(3 t . 
We find that in the limit [3 t — > the two-loop anomalous- 
dimension matrices reduces to 
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where the two-loop expressions for the anomalous dimen- 
sions 7cus P , 7 9 , and -f® can be found in [17], and 
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(12) 



arises from the threshold expansion of the two-loop coef- 
ficient of the velocity-dependent cusp anomalous dimen- 
sion 7 C usp(/3,a iS )- 

We stress that, as a consequence of the Coulomb singu- 
larities, the three-parton correlation term does not vanish 
near threshold. Instead, it gives rise to scattering-angle 
dependent, off-diagonal contribution in (11). The off- 
diagonal terms were omitted in two recent papers [22, 23], 
where threshold resummation for top-quark pair produc- 
tion was studied at NNLL order. We leave it to future 
work to explore if and how the results obtained by these 
authors need to be modified in light of our findings. 



V. CONCLUSIONS 

The IR divergences of scattering amplitudes in non- 
abelian gauge theories can be absorbed into a multi- 
plicative renormalization factor, whose form is deter- 
mined by an anomalous-dimension matrix in color space. 
At two-loop order this anomalous-dimension matrix con- 
tains pieces related to color and momentum correlations 
between three partons, as long as at least two of them 



are massive. This information is encoded in two univer- 
sal functions: F\, describing correlations between three 
massive partons, and / 2 , describing correlations between 
two massive and one massless parton. In this Letter we 
have calculated these functions at two-loop order. Using 
the exact analytic expressions, we studied the properties 
of the three-parton correlations in the small-mass and 
threshold limits. We found that the functions F\ and ji 
vanish as (mjmj / sjj) 2 in the small-mass limit, in accor- 
dance with existing factorization theoreoms for massive 
scattering amplitudes [14, 15]. On the other hand, and 
contrary to naive expectations, the two functions do not 
vanish in the threshold limit, because Coulomb singulari- 
ties compensate a zero resulting from the anti-symmetry 
under exchange of two velocity vectors. This fact has 
been overlooked in the recent papers [22, 23], where the 
three-parton correlations were neglected near threshold. 

Our results allow for the calculation of the IR poles 
in an arbitrary on-shell, n-particle scattering amplitude 
to two-loop order, where any number of the n partons 
can be massive. As an application, we have derived the 
anomalous-dimension matrix for top-quark pair produc- 
tion in the qq — > tt channel. We will explore in future 
work to what extent the new off-diagonal entries, arising 
from three-parton correlation terms, affect the numerical 
results for the threshold-resummed tt production cross 
sections at the Tevatron and LHC. 

This Letter completes the study of IR divergences of 
two-loop scattering amplitudes with an arbitrary num- 
ber of massive and massless external particles, and in 
arbitrary non-abelian (or abelian) gauge theories with 
massless gauge bosons. Details of our calculations will 
be presented in a forthcoming article. 
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